Inclusion of non-idealities in the continuous photodetection model 
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Some non-ideal effects as non-unit quantum efficiency, dark counts, dead time and cavity losses 
that occur in experiments are incorporated within the continuous photodetection model by using the 
analytical quantum trajectories approach. We show that in standard photocounting experiments the 
validity of the model can be verified, and the formal expression for the quantum jump superoperator 
can also be checked. 
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I. INTRODUCTION 

The continuous photodetection model (CPM) was pro- 
posed in the early 1980's in order to treat quantum optics 
situations in which a weak electromagnetic field enclosed 
in a cavity is continuously measured through the pho- 
tocounting approach The theory has received con- 
siderable attention in the following years due to its new 
microscopic interpretation of the photodetection process 
0, S HJPJHL relation to the quantum trajectories ap- 
proach [a 11, [13, [H [il and several proposals for ap- 
plications. Among them we find studies of photocounts 
statistics in diverse systems fisl. [ll . [isl ll ^ . . quantum 
non-demolition measurements lla. Il9l. 12011 . implementa- 
tion of measurement schemes 21 1221231. quantum state 
preparation f^, 25 , [1^, [13, HI] , quantum control via pho- 
todetection [23, 3C|, and quantum computation [3l| . 

CPM is extensively discussed in the literature [1, 
[3^ . [ssl . [3^ . so we shall mention only its main proper- 
ties. The model, also referred theory, describes 
the field state evolution during the photodetection pro- 
cess in a closed cavity and is formulated in terms of two 
fundamental operations, assumed to represent the only 
events taking place at each infinitesimal time interval. 
(1) The one-count operation, represented by the Quan- 
tum Jump Superoperator (QJS), describes the detector's 
action on the field upon a single count, and the trace cal- 
culation over the QJS gives the probability per unit time 
for occurrence of a detection. (2) The no-count opera- 
tion describes the field non-unitary evolution in absence 
of counts. 

If one sets the formal expressions for these operations, 
all possible outcomes of a photocounting experiment can 
be predicted. For instance, the photocounts [H, [2, Q and 
the waiting time [s^, [s^, [13, [11] statistics are among 
the most common quantities to be studied both theoret- 
ically and experimentally. Moreover, CPM conferred a 
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new step in photodetection theories by allowing to deter- 
mine the field state after an arbitrary sequence of mea- 
surements, thus creating the possibility of controlling the 
field properties in real time experiments [l6l . [itI . [soj . 

Actually, the QJS is the main formal ingredient within 
the theory, since it also dictates the form of the no-count 
superoperator Two different models for the QJS were 
proposed ad hoc. The first one was proposed by Srinivas 
and Davies the SD-model, as 



j p — Xapa 



(1) 



where p is the field density matrix, a and are the 
usual bosonic ladder operators and A is roughly the de- 
tector's ideal counting rate [l|, [11]. From the very be- 
ginning the authors [1| denounced the presence of some 
inconsistences when the QJS ([1]) is employed for describ- 
ing a real photodetection process, this point was also 
appointed in [3911. Nevertheless, this QJS is widely used 
in the literaturell, i, MM, [ll, M, M, M, 

The other proposal [3J, [40] assumes for the QJS an 
expression written in terms of the ladder operators = 
= E[_ (also known as exponential 



(a) a + 1) a and 



phase operators [4l|, 



Jp = XE^pE^ 



(2) 



In [3^ we called E-model such a choice, to differentiate 
from the SD QJS ([T]). Besides eliminating the incon- 
sistencies within the SD-model, the use of the E-model 
leads to different qualitative and quantitative predictions 
for several observable quantities. By an analysis of a mi- 
croscopic model for the detector, it was recently shown 
that the QJS's ((T|) and ^ are particular cases of a general 
time-dependent transition superoperator, each one occur- 
ring in a particular regime of the detector experimental 
parameters [H, [13]. Moreover, it was pointed out that 
by manipulating certain detector's accessible parameters 
one could engineer the form of the QJS, thus changing 
the dynamics of the photodetection, as well as the field 
state after a sequence of measurements. 

A way to check the validity of the CPM and to decide 
which QJS better describes the phenomenon in practice 
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can be accomplished through photocount experiments in 
a high finesse cavity by comparing the results to the the- 
oretical predictions. However, real detectors and cavi- 
ties are far from ideal. So our first goal is to include 
into the CPM the effects of non-ideality, such as quan- 
tum efficiency (QE), dark counts, detector's dead-time 
and cavity damping. Our second goal is to call attention 
to the fact that standard photodetcction measurements 
could verify which of the QJS models actually prevails 
exp er iment ally. 

The plan of the paper is as follows. In Sec. |TT1 we 
present a simple model, which enables us to include the 
effects of non-ideality - QE and dark counts - into the 
CPM using the quantum trajectories approach. Then 
we calculate the main quantities characterizing the pho- 
todetcction process - the photocounting and waiting time 
distributions. In subsection lll Al we do this using the QJS 
([I]), and in subsection lll Bl we repeat the same procedure 
for E-model. In Sec. IIIIl we analyze the behavior of the 
lower moments of the above distributions in realistic sit- 
uations and point out how one could decide about a QJS 
from experimental data. Sec. IIVI contains conclusions. 
In the appendix [X] we treat the effects of dead-time and 
cavity damping: we show that 1) cavity losses are not sig- 
nificant compared to non-unit QE effect, 2) the dead-time 
effect leads to mathematical inconsistences in SD-model, 
yet it is free of them in E-model, being however quite 
small compared to QE effect. Appendix [B] contains some 
mathematical details concerning evaluation of quantities 
of interest for different quantum states. 



II. MODELS OF NON-IDEAL 
PHOTODETECTORS 

A. SD-model 

We consider a free electromagnetic mono-modal field 
of frequency w, enclosed in an ideal cavity together with a 
photodetector (in the appendix[X]we show that the cavity 
damping is not crucial if the detector has non-unit QE). 
The unconditioned time evolution (UTE) of the field in 
the presence of the detector, i.e. the evolution when the 
detector is turned on but the outcomes of the measure- 
ments are disregarded (not reg istered), is described by 
the master equation [l^, [S^, [sj] 



p = —iuj (np — pn) — — iiip + pn — 2Ap 



(3) 



where Ap = apa' is a superoperator and ?i = a' a is 
the number operator. The first term stands for the free 
field evolution while the second describes the effect of 
the detector on the field due to their mutual interaction. 
The parameter A is the field-detector coupling constant, 
roughly equal to the ideal counting rate [4^, [43] . 

To describe photocounting with QE rj and finite dark 
counts rate Xd (d is the ratio between the dark counts 
rate and the ideal photon counting rate), we assume the 



following expression for the QJS (c.f. the expression re- 
sulting from the microscopic model in [4?}) 



(4) 



Jp = X(j]A + d] p 



It describes the action of the detector on the field upon 
a photodetcction, and its trace gives the probability per 
unit time of the click. Actually, the microscopic model 
[46l | suggests that Q has a diagonal form in the Fock 
basis, but this will not be important here, since we 
shall be interested only in diagonal elements. The first 
term within the parenthesis describes the absorption of 
a photon from the field with probability per unit time 



Tr 



rjXAp 



= rjXn, where n is the field mean photon 

number - this means that the detector 'sees' all the pho- 
tons. The second term describes the occurrence of a dark 
count with field-independent probability density Ad, and 
this event by itself does not modify the field state (the 
field state after a single dark count is Xpd/Tr {Xpd) = p). 
However, when both terms are present, the field state 
upon a detector's click becomes a mixture of both out- 
comes. 

From the quantum trajectories approach and CPM 
[l|, 0, 01 , all the quantities related to photodetcction can 
be calculated provided the complementary no-count su- 
peroperator St is known (St describes the action of the 
detector on the field during the time interval t without 
registered counts). Acting St on the initial field state 
Pq, the no-count state pg = Stp^ obeys Eq. ([3]) when one 
subtracts the term (g]) on the RHS (see (EH)- Moreover, 
as we are interested in calculating probabilities, we shall 
disregard phase factors cxjp{±iujnt) , since they are can- 
celed in any trace evaluation. So the evolution equation 
of pg is 

pg ^ -^{hpg + pgn) + XqApg- Xdpg, q = l-r]. (5) 

Setting the transformation 

p^ = e-''^*Utp„ Utp = e-^*"/2pe-^*"/2 (6) 
in Eq. ([5]) we obtain a simple equation for p^ 

Pi = Xqe-^'Ap^, (7) 

whose solution is 



Pi 



E 

/=0 



^Vo(«0' =cxp(#ti)Po- (8) 



where 



1 - e 



-At 



Thus the no-count superoperator is 



StPo 



-dXt 



(9) 



(10) 
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The field UTE supcroperator Ti, defined as the solution 
to Eq. ([3]), is naturally given by setting d = 77 = in 
Eqs. dS]) and i.e. 



ft = Utie'^^^Po). 



(11) 



We introduced in Eq. ^ a compact notation for the in- 
finite sum in terms of the exponential supcroperator. We 
can deal with such superoperators as they were common 
operators, provided we use the 'commutation relations' 



AUt = e-^^UtA 



ey^Ut^UtCxpiye~^'A), (12) 



obtained by expanding the superoperators in series. 

Now we can calculate the m-counts supcroperator 
Nt{m), that describes the field state after m registered 
counts (whatever real or dark ones) in the time interval 
(0, t), and whose trace gives the probability for this event. 
It reads 



Nt{m)p- 



hp, 



(13) 



where the integrals are evaluated over all the time inter- 
vals between the counts 



= I dUa / dt„,-i ■ ■■ dh (14) 
Jo Jo 



and the conditioned density operator is 

hp=St^t^JSt^-t^_^J---JStiP. (15) 

Expanding the QJS ([4]) in Eq. (fTS)) in terms of rjA and 
d, one obtains a finite sum whose first term, proportional 
to describes the detection of m photons: 

= {\r,r St-t^A- ■ ■ ASt, (16) 
= (A77)™e-^(*^+*^+-+*")^ti™. 

After integrating we obtain the first term in ([13]), 
describing the field state after the loss by absorption of 
m photons, 



fit{m) — 



(17) 



Calculating in a similar way the contribution of the terms 
with higher powers in d we arrive at the formula 



iVt(m) = ^ 



fc=0 



{dxty 



ht{m — k) = St 



{dXt + r/0ji)" 



(18) 

One can easily verify that the m-counts superoperators 
(dH]) satisfy identically the fundamental relation P, 0] 



Nt{m)=ft. 



(19) 



The factorial moments of the photocounts distribution 
are easily evaluated as 

00 

TO • • • (m — l)^ ~ TO • • • (m — /)Tr[7Vt(m)/9] 

= Tr[5t(dAt + r]0ti)'+^ exp(dAt -f ?y(^ti)p] 

= Tr[J7t(dAi ^ 7?0ti)'+ie'^*^p]. (20) 

Thus we need to calculate the expression 
$fe(6,2;) = Tr[j76e"^iV 



E 



n\l\ 



{n-k)\ 



(x 4- e 



-\b 



where p„ ~ {n\p\n). Evaluating 

00 



ri=0 



(n - fc)! 



(21) 



(22) 



(see Eq. © for the expression of (pf ) we obtain general 
expressions for the lower factorial moments 



TOt = dXt + rjucji^ 



(23) 



to(to - 1)( = {dXt) 4- 2r]ndXt(j)t + {r](j)tyn{n - 1), (24) 



where h and n{n — 1) are the factorial moments of the 
initial density operator. 

Another measurable quantity we consider here is the 
waiting time distribution. It describes the probability 
density for registering two consecutive clicks separated 
by the time interval r, under the condition that the first 
one occurred at time t. Its non-normalized form is 



Wtir) = Tr[jSrJTtp 
and the mean waiting time is 



--Af' 



drWt (r) r, TV : 



drWt (t) 



(25) 



(26) 



where T is the time interval during which one evaluates 
the averaging in experiments. As will be shown in section 
mil T is an important parameter due to the presence of 
dark counts. After straightforward manipulations, using 
the 'commutation relations' (|12p. we obtain 



Wt(r) 



+ride-^\l 



dH^] 



where 



<l>r = <i>4* + 1 - e-^* {v+il- 



(27) 
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In the appendix [X] we consider the dead-time effect 
and show that it cannot be consistently incorporated 
into SD-model, because the QJS (HI is an unbounded 
superoperator and the resulting counting probability is 
non-normalizable. This is just one more mathematical 
inconsistency of the SD-model. In the appendix |B] 
we evaluate the expression (PT|) for three kinds of states: 
coherent, number and thermal. 



and we note again that all the composite superopera- 
tors, such as exp(?/e), are understood as power expan- 
sions. Now, multiplying (pij) by Aq on both sides and 
using the solution (|33p we get an equation for AoPj^Aq 



di 

with solution 



(AoPiAo) = XqAo (e-^*(i-«-')epo) Aq (34) 



B. E-model 

We now repeat the same procedures for E-model in 
which the QJS is 



Jp = X {'qe + d) p. 



(28) 



where ep = E-pE^. The probability per unit time for 
detecting a photon is 77A(1 — Po), where po = (0|p|0), so 
the detector 'sees' whether there is any photon in the 
cavity. In principle, the parameter A is different from 
the one in SD-model, but in the context of this paper it 
will be always clear which one we are dealing with. The 
field UTE is described by an equation similar to Eq. ([3]), 
obtained by doing the substitution |a,a^} {£'_,i?+} 
in the non- unitary evolution (second term on the RHS). 
So the no-count state pg obeys the equation 

[similar to Eq. (O] where A = E+E^ = 1 - Aq, Aq = 
|0)(0|. Setting the transformation 



AoPiAo = Ao 



1 - qeR t 
1 — qe 



-Po 



Ao, i?t = e-^*(i^«^). 



(35) 

Thus the diagonal form of the no-count superoperator, 
which we write just in terms of the projector Aq and the 
unit operator, is 



S, 



tPo 



-dXt 



Rt + Ao- ^Ao 

I — qe 



(36) 



where we use the notation (AoQAq)/? = Ao{Qp)Ao. 

Repeating the steps - (UHl), we obtain first the 
conditioned density operator 



p^e 



-dXt 



Rt + Ao^^f^Ao 

1 — qe 



{XfieTp. (37) 



After evaluating the time integrals as in pT)l we get 



fitirn) 



1 — qe J m\ 



+Ao 



{hfe 



„m — 1 



dxRx- 



l~qe Jo (m- 1)! 



rAo 



(30) 



in Eq. (^9)) and using the property exp(aA) ~ Aq -I- e"A, 
we obtain the differential equation for pi 

Pi = Xqe~^* (Ao + e^*/2A) E^ (ApjA) (31) 
x^+(Ao + e^*/2A). 

We solve this equation by projecting it onto orthogo- 
nal subspaces spanned by projectors {A, Aq}. Moreover, 
since at the end we shall be interested only in calculat- 
ing probabilities, we consider only the diagonal part in 
Fock basis for quantities of interest, thus disregarding the 
terms whose trace is null, such as ApA^. 

Multiplying Eq. (|3ip by A on both sides we obtain 



-{Ap,A) =XqAE^[Ap,A)E+A, (32) 



for m > and ht{0) = St- Finally, analogously to the 
expression (|18p. we obtain the m-counts superoperator 



Ntim) 



fAo- ^Ao 



1 — qe ml 



(38) 



-Ao 



Xr]e 
1 - 9^ 7o 



dxRx 



[dXt + rjiXx 
(to - 1)! 



-Ao 



where the last term is zero for m = 0. One can easily 
verify that the superoperator Nt{m)^ Eq. ([55]) . satisfies 
relation p9)) . 

After lengthy however straightforward calculations we 
obtain the following expressions for the initial factorial 
photocounts moments 



rrit = dXt + 7]n (1 — ^i) . 



(39) 



whose solution is 



Ap,A = A (e'^*Vo) A 



(33) 



to(to - l)j = [dXt) + 2i]ndXt (1 - Si) 



(40) 



+ -q-^ n{n - 1) (1 - fi) - 2nXtE2 
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where 



-Tr 



i?0 = Rtiq = 1), (41) 



n = 



n{n — 1) 



=Tr 



1-e 



(42) 



Using Eq. (pS)) . the waiting time distribution density is 
found to be 



TiiiJRr + AdAo^— ^Ao)Ji??p] I . (43) 



In the appendix [Xl we show that the dead-time effect 
can be incorporated into E-model, however its effect is 
quite small compared to the non-unit QE effect, so we 
disregard it in this paper. In the appendix IB] we obtain 
formulas for Eqs. (gT]), ([42]) and (gSl) in terms of p„ 
and evaluate them for the coherent, number and thermal 
states. 




200 300 



500 



FIG. 1: Mean photocounts number fht in E-model for coher- 
ent, number and thermal states (indicated in the figure, the 
lower curves are labeled analogously) as function of time for 
two values of the initial photon number: the lower curves cor- 
respond to n = 50 and the upper - to n = 100. In the inset 
we plot fht for the SD-model, which is independent from field 
state. 



III. VERIFYING CPM 

Basing ourselves on published experimental data (48l | 
we chose the following numerical values for the model 
parameters: 77 = 0.6 for the QE and d — 5 ■ I0~^ for 
the dark counts rate (normalized by the ideal counting 
rate). We do not attribute any fixed value to A since 
our analysis will be given in terms of the dimensionless 
Xt. Many photodetection quantities in different contexts 
were reported in, e.g., [l|, H, [13, [13, [H, [H, [H, [11] , so here 
we shall consider few of them that could check the valid- 
ity of either, the SD- or the E- model in photocounting 
experiments. 

First we analyze the counting statistics. In figure [1] we 
plot fht as function of Xt for both models for two val- 
ues of the initial mean photon number, n = 50 and 100. 
Initially, 7714 increases steeply due to photons absorption, 
and after some time the growths turns linear with much 
smaller slope due to the dark counts. We call the time 
interval during which the photons are absorbed (repre- 
senting the duration of the steep increase in the number 
of counts) the effective counting time tE- In the E-model 
t^; is proportional to the initial average photon number, 
contrary to the SD-model [as seen from the figure [1] and 
formulas ((23)) and (|39p ]. So the experimental analysis of 
the dependence of < ^ on n seems to us a feasible way for 
verifying which model could hold in practice, because, 
according to the SD-model, ts does not depend on n. 
Moreover, one could also check the validity of each model 
by verifying whether fht depends on the initial field state: 
in the SD-model it is independent of the field state, while 
in the E-model fht is quite sensible to it: in figure [T] one 



sees a notable difference between thermal and coherent 
states, although not so much between number and coher- 
ent states. This can be explained by a great difference 
in the values of Mandel's Q-factor [i^] characterizing the 
statistics of photons in the initial state: it equals —1 and 
for number and coherent states, respectively, whereas 
it is very big {Qth = fi) for the thermal states with big 
mean numbers of photons. 

Now we analyze the normalized second factorial mo- 
ment 



Kt 



(44) 



for the same initial states with mean photon number 
n = 50. For the number and thermal states Kt as func- 
tion of Xt is shown in figure [21 and for the coherent state 
we get Kt = 1, so it is not plotted. In the asymptotic 
time limit and for non-zero dark counts rate, the same 
value Kao ^ 1 holds for both models, however the tran- 
sient is model dependent. In the SD-model without con- 
sidering dark counts Kt is time- independent, writing as 
K ~ fi{n — l)/«^ [h and n{fi— 1) correspond to the ini- 
tial field state), nevertheless it depends on the initial field 
state: K = 2 for the thermal state and K = \ — \/h iox 
the number state. By including the dark counts in the 
analysis this constant behavior is slowly modified as time 
goes on, see figure [2] 

In the E-model in the absence of dark counts Kt starts 
at the value 



limi^:, 



Tr 



1 - Po - Pi 



(45) 
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FIG. 2: Normalized second factorial moment Kt, Eq. (|44|l . 
for SD- and E- models (as indicated in the graph with abbre- 
viations) for the number state (and the thermal state in the 
inset) for n = 50. For the coherent state one has Kt = 1 a,t 
all times for both models. 



FIG. 3: Mean waiting time ft as function of Ncav for the 
number (N) and thermal (T) states for SD- and E- models. 
While there are photons in the cavity ft is constant for the 
E-model, but increases with time for the SD-model. In the 
inset we plot Ncav as function of \t for these states (in the 
SD-model Ncav is state independent). 



which is exactly 1 for the number state and very close to 
1 for the thermal state with the chosen values of n. With 
the course of time, Kt attains the same values as for 
the SD-model (for respective initial field states) when all 
the photons have been counted. By taking in account the 
dark counts effect such a behavior is slightly modified, yet 
it is quite different from the behavior in the SD-model, 
as shown in the figure [2l This is another possible manner 
for verifying the applicability of SD- or E- models. 

We now turn our attention to the waiting time analy- 
sis. It is important to define the time interval in which wc 
do the average: if one has non-zero dark counting rate, 
then by performing the average over a very large time in- 
terval, we shall always get for the mean waiting time the 
value f ~ (Ad) ^, which is nothing but the mean time 
interval between consecutive dark counts. Since experi- 
mentally the average is done over finite time intervals, 
we shall proceed in the same way: the mean waiting 
time for initial times, when the photon number is sig- 
nificative, is roughly {t]X)~^ (because rjX is the effective 
counting rate), so wc shall take the average over a time 
interval ly — 10 (ryA)~^. This means that if one docs not 
detect consecutive counts within the time such a mea- 
surement will not contribute to the average. In an ideal 
case this procedure is not necessary because the prob- 
ability for registering consecutive clicks separated by a 
large time interval is zero. 

In figure [3] we plot the mean waiting time for the SD- 
and E- models, for the number and thermal initial states 
(for the coherent state we obtain a curve almost identical 
to the one for a number state) with n = 100 as function 
of the mean photon number in the cavity at the moment 



of the first click. 



CAV 



= Tr 



for SD-model 



n^i for E-modcl. 



(46) 



(For completeness, in the inset of figure [3] we plot Nqav 
as function of Xt for both models.) For the E-model, 
we see that when Nqav becomes less than 1, the wait- 
ing time starts to increase dramatically due to the dom- 
inance of dark counts, which are much more rare events 
than absorption of photons. This is a drastic difference 
from the ideal case, in which no counts occur after all 
the photons having been absorbed, so the mean waiting 
time saturates at the inverse value of the counting rate, 
as shown in [3^. Moreover, from figure [3] one verifies that 
as long as there are photons in the cavity the mean wait- 
ing time is nearly time-independent within the E-model 
(and truly independent in the ideal case (soj). and do 
increase substantially in time for SD-model. This is an- 
other notable qualitative difference we suppose one could 
verify experimentally. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have generalized the continuous pho- 
todetection model through a careful quantum treatment 
of non-ideal effects that are ubiquitous in experiments. 
We derived general expressions for the fundamental oper- 
ations in the presence of non-unit quantum efficiency and 
dark counts, and calculated explicitly the photocounts 
and the waiting time probability distributions for initial 
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coherent, number and thermal field states. By calculat- 
ing the first and second factorial moments of the pho- 
tocounts and the mean waiting time, we showed that in 
standard photodetection experiments one could check the 
applicability of the QJS of SD- or E- models. Namely, 
we indicated three different ways for revealing the actual 
QJS: (1) quantitatively, one should study the time de- 
pendence of the normalized second factorial photocounts 
moment. Qualitatively, we showed that the models can 
be also distinguished by measuring: (2) whether the effec- 
tive detection time depends on the initial average photon 
number in the cavity and (3) whether the mean waiting 
time is modified as time goes on. To that end we have 
considered three different kinds of field in the cavity: the 
number, coherent and thermal states; each one on its 
own permitted to do comparisons between the two stud- 
ied QJS's. Results with other kinds of fields could also be 
presented here, as for instance the binomial state or the 
so-called squeezed state, however, no new physics related 
to the goals of the paper appears. A last remark, if the 
experimental data would depart significantly from the 
theoretical prediction one should reconsider both models 
and try to look for alternative mechanisms to reproduce 
the outcomes. 

In conclusion, we believe that our theoretical treatment 
could provide clues for an experimental verification of 
the CPM, contributing with valuable insights about the 
quantum nature of the photodetection in cavities, as well 
as giving rise to the possibility of field state manipulation 
through detector post-action on the field. 
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APPENDIX A: CAVITY DAMPING AND 
DEAD-TIME 



p=l 



q = 1 — 77 + c = p — 77. 



(A3) 



The value of c should be at the most of order of 10^^ in 
order to make viable the CPM. In this case we see that 
if one takes into account the QE drawback, the cavity 
damping does not modify substantially the resulting ex- 
pressions. Therefore wc disregard its effect in this paper. 

The dead-time effect means that immediately after a 
click the detector is unable to register another count 
within a quite small time interval x. Ax <C 1. In our 
framework we can describe this effect as the occurrence of 
the UTE during the time x immediately after the count, 
so the conditioned density operator h p3|) becomes 



hp = St-t„ 
= St-t„ 



QJSt^.t^_,QJ ■■■QJSt.p, (A4) 



where the dead-time superoperator, under condition 
Ax ^ 1, is found to be 



Q = S-xTx ~ cxp(.Tj) 



(A5) 



for both SD- and E- models, with respective QJS's. 

In SD-modcl the resulting dead-time superoperator is 
unbounded, as well as J, so it can bring some mathemat- 
ical inconsistences. For example, the m-counts superop- 
erator with dead-time effect is found to be 



Nt{m) 



[dAi-He'*^^z/(p(/)J]' 



where 



(A6) 



(A7) 



one 



If one evaluates, for instance, Tr 

will find a divergent result because z increases much 
faster than the decreasing terms. 

In the E-model J is a bounded superoperator, so the 
dead-time corrections will be of order 77 Ax <C 1, much 
less relevant than the non-unit QE drawback. 



APPENDIX B: EVALUATION OF TRACES 



First we include the effect of cavity damping in our 
treatment. In quantum optics experiments the back- 
ground photons number is negligible, so we can model 
the cavity as a thermal reservoir with zero mean excita- 
tions number, described by the standard master equation 
0. Then the UTE equation in SD-model should be 



-iuj {hp — pn) — — ( hp + ph — 2Ap 
Ac 



(jip + ph — 2Ap^ , 



(Al) 



where Ac is the cavity damping rate. From it, following 
the steps of sec. |TT]we obtain the no-count superoperator 



(e'^^'^Po) , 0t = - (1 - e-^P*) , (A2) 



In this appendix we derive general expressions for both 
SD- and E- models and evaluate them for a general initial 
density operator p = X) Pnl'^) ('^1 (non-diagonal elements 
do not contribute to the trace in the expressions below). 
We shall analyze three particular field states: coherent 
state, 



Pn 



number state, 

Pn = ^■ri,n with integer fi, 
and thermal state, 

p„ = (l-a)Q!", a = n/(n+l), n(n - 1) = 2n^. 



In the SD-model. formula 



results in: 
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• coherent state 

$}f = n'' exp [-7;n0^e-^*] (Bl) 

• number state 



71 — k 



thermal state 



fc!(l - a)a'' 



[l-a(l-#^e-^*)] 



(B2) 



(B3) 



• For the coherent state we use the formula 



E 

k=0 



k\{k + n)\ 



where Ik (x) is the modified Bessel function [50| , to 
obtain 



^ , ^ . (n+fc)/2 

^(n+l)(-) /„+,(2VnAi), 



(B9) 



The formula (gS]) yields Ncav = ne~-^^ for all the 
states. 

In the E-model we need to expand the superoperators 
as series of e and evaluate the sums. For Eqs. ([1T|) and 
(l42l) we obtain 



E 



(At)' 



n — ' m 



-At °o 



I Pti+l+m+k 

(Xt) 



J2 {n + l)^-^Pn+m+k (B4) 



2e 



n.m— 



n(n — 1) , „ 

„-At °° 



" ml ^"+'+"1+1 



n(n — 1) 



E "(^-1) 



m 



(B5) 



Regarding the evaluation of the mean waiting time (|43|) . 
one needs to evaluate the expressions 



Tr 



An 



1 — qe 



(B6) 



— At — n - rj/2 

n=^==Y,n{n-l)(^\ In{2VfJd), (BIO) 



ji=0 ^ ^ 

(Bll) 

The above series can be transformed in a finite in- 
tegral using 



°° ptz/2 fZ 



valid for Re(i^) > 0. 

For the number state, using X]fc=o ^'^ 1^^- = e'^r(n-|- 
l,a;)/n!, where r(a,a;) = /^t"~^e~*dt is the in- 
complete complementary Gamma function [soj , we 
obtain 



r(?l - fc + 2, At) - Air(7i - fc + 1, AO 

^k = rrr T-r. , (B12) 

nyn — k)\ 



where 



^k{q.P)^ E 9"^Pn+;+fc- (B8) 



.,1=0 



(B7) 



For the thermal state we can evaluate the expres- 
sions obtained in the section Hi Bl directly using the 
'eigenstate' relation ep = ap and Tr[AopAo] = p^. 



T{n + 1, Xt) - 2Atr(n, Xt) + {Xt)'^r{n - 1, Xt) 



i{n- l)(n- 2)! 



(B13) 



*.(.,/3) = ,^-'^e^^/.£S^-h±l^. (B14) 

[n — k)\ 
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